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CN . Abstract 

Variation principle has been developed to calculate many-particle effects in 
crystals. Within the framework of quasi-particle concept the variation princi- 
00 \ pie has been used to find one-electron states with taking into account of effects 

due to non-locality of electronic density functional in electromagnetic fields. A 
secondary quantized density matrix was used to find the Green function of a 
quasiparticle and changes of its effective mass due to correlated motion of inter- 
\ acting electrons. 

§ : 1 Introduction 

o. 

Because of high dimensionality the complexity of many-electron problems all problems of a 
solid state physics of this type are much higher than of one-electron problem. Results of 
a solution of the one-electron problem can be utilized in many-electron problems in solid 
state physics if one supposes that motion of an electron happens in a self-consistent one- 



particle potential V(r). The self-consistent potential is yielded by the solution of the Poisson 
equation, for which the density of one-electron states is a self-consistent solution of the one- 
particle Hartree - Fock equations |T] - [3]. These equations have been written originally 
to calculate one-particle states of many-electron atom. Equations of the Hartree - Fock 
type for band calculations have been proposed by Kohn and Sham jS]. It is practically 
impossible to look for their solution without additional assumptions if many-particle effects 
are circumscribed by distributed in space electron density functional. One of them is that 
one-particle solutions describe a crystal as a set of interacting quasiparticle excitations 
In a given paper we offer an approach which allows to determine one-electron states with 
taking into account of the effects owing to a non-locality of electron density functional in 
electromagnetic fields and which is based on usage of the variation principle. 



2 A variational principle for band calculations 



Let us write a one-electron Hartree - Fock Hamiltonian for a system of N electrons, N 
moving in a field originated by atoms nucleus' of a crystal 



oo, 



:Pi + U(n) + V sc (n, a t ) - t x {r u a t 



"tpnik^i 



(e n (0) + e„(/c i ))V'„(A; i r i ) 
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with a unit choice h = 1, m = 1, where pf/2 = — |A(f\) is a kinetic energy in system of 
atomic units, A(r*j) is the Laplacian operator written in a given point with a radius- vector r*j 
in which i-ih electron having a spin <7j is situated; £7(rv) is potential energy of i-th electron in 
nucleus field of the crystal, V sc and S x are operators of Coulomb and exchange interactions, 
respectively 



N 



m=l 



N 



m=l 



(2) 

S x (f i ,cr i )V'n(fci^) = ^21 Vmihr'^v^ri- fJ|)^„(A*rJ) dr^^n); 

(3) 



rj = {r*j,o"i}, ipnikiTi) is a wave function including spin and coordinate parts, u(|rj — r*|) is 
potential energy of electron interaction. The physical sense of the operators (0 El) becomes 
obvious if one rewrites them in terms of spin-zero electronic density p(r, r") and assumes 
that the interaction v is the Coulomb one: 



N-l 



(jV-l)/2 



P(r, O = j (^m(^ • f ' • r , -a) + ^(fc • f, -a)i[) m (k ■ 7, a) 

N-l)/2 

ip* m (k-r)ilJm(k-r'), v = e 2 /\r-7\. (4) 



m=l 



It follows from here that the operator V sc gives electrostatic interaction of one electron with 
an electron density created by residual N — l electrons, with an electrostatic self-action (s.a.): 

(iV-l)/2 

V sc {f u a^nihn) = / ■ ?M\ri ~ KD^mih • rj) d^ n {h ■ n) + s.a. 

a m=l 

d ^^^ Mhrd+s.z. (5) 



in - r! 1 



Analogously we get that the operator gives a quantum exchange with an exchange self- 
action (s.a.): 

N ^ f f 

t x {f u o-j)V'n(^n) = 22 I drjdajijj* m (ki ■ fj, (Tj)v(\fi - f^-|)^ n (A* • r}, <Tj)if) m (ki ■ f u a { ) 

1 iV_1 /" f 



m=l 



m=l 



+i)* m {k ■ fi, -aj)ij m (k ■ fj, (Jj)S{crj - ov)J f (|f; - r^|)^ n (fcj • r}, <Tj) + s.a. 



e 2 p(r j ,r i ) 

I TV - 7\j 



d?j ~ Z v J, ^Y ipnikrj) +s.a. (6) 



Since operators V rsc and E x in expression (JIJ are subtracted one from another the self-acting 
terms vanish. 
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A quantity e n (ki) entering in expansion E n {ki) = /(e n (0)) + e n (ki) being n-th eigenvalue 
E n (ki) of the Hamiltonian for N particles system in Hartree -Fock approximation is the 
energy of n-th band. A quantity e„(0) determines a reference point ExtiE n (ki) = /(e n (0)) of 
n-th bands, where / is an unknown function. Within the framework of quasiparticle concept 
the physical sense of e n (kj) is the energy of quasiparticle excitation. 

Let us take a solution of one-electron problem for an atomic area in a cell potentials 
approximation as e n (0) and a basic set on which one constructs an expansion of a trial 
function ipnikiTi) of n-th band. Then one can realize a variational principle as the following 
expression: 

5e n [^ n ] = 0, (7) 

which is determined not strictly as it is used to find excited states for which, as a rule, the 
basic set of functions is unknown 

Let us show, that the variational principle for excited states in the form (J7J) can be made 
strictly determined one. Since, by definition, the energy of a nonexcited state coincides with 
an extremum of a functional E n [if) n ] we have the following equality of variations: 

5eM = 5E n $ n ] = 0. (8) 

The ambiguity of a variation (jSJ) consists only in an arbitrariness of a position of a reference 
point for an energy band as the variation procedure for an one-electron problem is strictly 
determined. If a symmetry of atomic areas is definite and a cell partition of the crystalline 
space is unambiguous then, in principle, this arbitrariness is removed easy by utilizing the 
quasi-particles concept, according to which 

£ n {ki) = e n (-ki). (9) 

It follows from the expressions (jHJ) and 0, that the additional variation of a reference point 
of a band gives a coincidence of this point with a centre of the energy band for a crystal and 
allows to transform the variational principle (jHJ) finding noninteracting quasiparticles states 
to strictly determined expression 

5E n $ n ] = 0, 8f$ n ] = 0. (10) 

Let us find the reference point ExtrE n (ki) of energy band using a method of density matrix 
functional. 

3 Secondary quantized reduced density matrix 

As is known jH] , equations in quantum mechanics can be written not for a wave function, but 
for a density matrix p. The operator p is a projective operator for pure states and can be 
presented in terms of Dirac ket(bra)-vectors as p = If the operator p is known then 

we can find, by definition, energy E of system of N particles, described by a Hamiltonian 
H(fi, . . . , fjv) and a wave function as 



E = Sp pH. 



(11) 
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Let us introduce n-dimensional density matrix p n by 



Pn( r \i ■ ■ ■ i r n'i r l' • • • ' r n) 



(N-n)\ 



dr 



n+l 



dr 



N 



X (^l> ■ ■ ■ y n ^n+l, ■ ■ - ,7^1-0X^1^1, • • .,r n ,f n+1 , . . .,f N ), 



(12) 



which, by implication, is a reduced coordinate distribution n-particle function [Oj. By defi- 
nition, it has a normalization 



SpPn = TTT 77 U = 1, . . . , N. 

(i\ — n)\ 



(13) 



Let us examine the system consisted of N with a pairwise interaction 

JV 
i>j=l 

Then, since projective operators possess the following properties: p 2 = p and p* = p the 
expression (jllj) is transformed to the form: 



N 



N 



N 



e = Spp(h + v (\n - fj\) J = Sp p^2 Kn) + Sp p 2 Y v ^ *v) = 

V i>j=i J i=i i>i=i 

JV TV 
8=1 i>j=l 

where h(fi) = |p| + Ufa). Since the non-perturbed hamiltonian ifo m Eq. ()14|) consists of 
independent one-particle summands and interaction of particles occurs pairwise, the operator 
of trace appears in Eq. (|14p over one variable r\ and two variables r 1; r 2 : 



= Sp M^i) / ^2 • • • dr N 
i=i 



N 



1 

X- 

2 



x{r[,f 2 , . . . ,r N \ip)(ip\fi, . . . ,f N ) +Sp u(|fi - r 2 | 

i>i=i 



rfrirff 3 . . . dr N |(fi,f^,r 3 , . . . , f/v|^)(^|fi, . . . ,r N )\' 



(15) 
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Using the equality: J k ■ kdO = 2\k\ 2 one can transform Eq. (|15J) and obtain: 

E = Sp h(fi)N / df 2 . . . dr^ 



x(r' 1 ,r 2 , . . . ,r N \il))(il)\r u ...,r N ) + ~ Sp v(\r x - r 2 \)N(N - 1) 



x / dfidf 3 ... df N {f 1 ,r 2 ,f 3 ,...,f N \-ip) x 
x I df 2 ... dr N {ip\r[,r2,r 3 , . . . ,rjv)(fi, • • • ,r N \i>)(i/>\ri, • • • ,r N ) 

+ J dr 2 . . . df N {r[,f 2 ,r 3 ,. . . ,r N \ip) x 
x / dndf 3 ... drNiipln,^,^, ..^rjvXfi, . . . ,r N \ip)(il)\ri, . . . ,r N ) 



(16) 

It is easy to see, that the first summand from the right-hand side of Eq. (JTfij) contains 
the reduced one-particle density matrix px as a multiplier. Therefore after some obvious 
transformations Eq. can be rewritten as: 

E = Sp h(n)p 1 (f^ 1 ,f 1 ) + ^Sp v(\n - f 2 \)N(N - 1) 
xj dr 3 ... dr N (^ 1 ,f^ 2 ,f 3 ,...,f N \^)(^\r 1 ,f 2 ,f 3 ,...,r N )x 

x J dr x dr 2 ... dr N {r x ,...,r N \i>)l^)\r x ,...,r N ) 
= Sp h(f 1 )p 1 (r[,f 1 ) + ^Sp v(\n- f 2 \)N{N - 1) 
x / df 3 ... df N (r[,r 2 ,f 3 , ... ,r N \i}>) (V>|fi, f 2 , f 3 , ... ,f N ). 



(17) 

Here one takes into account the normalization of wave function j dr\dr 2 . . . dr^ = 
1. 

If we observe that in the right side of the equation (|T7j) the second summand contains 
the reduced two-particle density matrix p 2 as a multiplier then it is possible to transform 
this equation to the expression which has the reduced density matrixes p\ and p 2 : 

E = Sp h(r 1 )p l (r' l ,r 1 ) + ^Sp v(\n - f 2 \)p 2 (f{, f^; ri, f 2 ) = 

= e {0) N + is P v(\n - f 2 \)p 2 (^y 2 ; n, f 2 ), (18) 

where is the one-electron state, the summand E HF = |Sp v{\f\ — r 2 \) p^r 1 ^ f^; fi, r 2 ) 
is excitation energy of the system under the interaction. Since in the Hartree - Fock ap- 
proximation the reduced density matrix p 2 is factorized, at first, the energy e^N equal to 
energy of N one-electron states, including kinetic energy of an electron, energy of an electron 
in a self-consistent scalar potential, and exchange energy, and, secondly, energy E HF of the 
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excitation yield, as it follows from the equations Q and (JIEJ), a contribution to the electronic 
energy E of a crystal. Further we shall show, that the Hartree - Fock approximation is an 
one-particle approximation in the sense that in this approximation the excitation energy 
E HF is represented as the energy of quasiparticle states. 

Now, we rewrite Eq. (P) in the representation of Dirac ket(bra)- vectors: 



h(k)\n; k) + Yl J S(k - k') dk'(\n;k)(m;k'\v(kk')\m;k') + 

m=l,mj^n 

+(|n; k')S nm )(m; k'\v(kk')(\m; k)5 mn )) 

N f 

— / \m;k)(m;k'\v(kk')\n;k')6(k — k')dk' = 

m=l J 

= \n;k){e n {0) + e n {k)) (19) 

where ki = {fcj,CTj}, h(k) is a momentum representation of the non-perturbed hamiltonian, 
v(kk') = J dfdr'\r)(k ■ f\v(\f— -P' | ) | A;' • r > )(r > \ is a momentum representation of the Coulomb 
interaction operator, 5(k — k') is the Dirac 5-function manifesting the presence of the law of 
conservation of momentum. 

Let us introduce projective operators p™£ 

p^ = \m;k')(n;k\ (20) 

and express Eq. (fT9~|) via these operators p™, . For this purpose, Eq. (fTT?|) is multiplied on the 
right by bra- vector (n; k\. Then, additional summating over n and integrating over dk one 
get the equation: 

N N r r 

I /(-/i i/\ -tit 

X 



[ dk h{k)\n- k){n; k\ + ^ / / 5(k - k') dk dk' 

n=l J n=l J J 

■k)J2(m;k'\v{k,k')\m]k')(n-k\+ j j 5{k - k') dk dk' x 

m=l J J 

N \ / N \ 

Y l n ! k')5 nm (m; k'\ v(kk') ^ l m ? k)S mn (n; k\ - 




\m=l / \n=l 

N N 



K k)(m; k'\v(kk') \ n ; k ')( n \ k \ 5 ( k> ~ k)d(-k)d(-k') 

N 

dk^2(n;k\n;k){e n {0) + e n {k)). 



m=l n=l 

N 



n=l 



(21) 



We see that the first and second terms on the left of the equation fj2*Tj) are traces of a 
matrix representation of operators ph and p*pv, and the third and fourth terms on the left 
of equation (J21j) are mutually cancelled. Hence, it means that using normability of function 
|n; k): f dk(n; k\n; k) = 1, we obtain the following equation: 

N 

Spph + $pp*pv = e n (0)N+ dk^2(n;k\n;k)e n (k). (22) 

J n=l 
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Using properties of the projective operators p™$: {p™$)* = and {p™$) 2 = p^§ one can 
transform Eq. f!22l) to the form: 

AT 

Spp(fr + v) = e n (0)N + / dfc ^2(n; k\n; k)e n {k) = e n (0)N + e. (23) 

n=l 

Let us elucidate a physical sense of introduced projective operators p™$. It follows from 
comparison (fTHj) and (J2HJ) that the energy of a quasi-particle e is on the right of Eq. (J23J) 
accurate to the constant e n (0)iV. It follows from here that the operator p™£ allows to 
calculate energy e of quasiparticle excitations. It means that the expression (l23|) is nothing 
else but a procedure of average on density matrix. Since the averaging with the help of the 
operator p™$ yields energy e of quasi-particle this operator is a secondary quantized density 
matrix. 

It follows from the comparison of right sides of the equations (JTJ), ([T%]). and that the 
reference point ExtrE n (ki) of energy band determines the solution e n (0) of the one-electron 
problem 

e n (0) = Extr E n {h)/N. (24) 

Thus, one has proved that the equation (JTJ) can be considered as an equation describing 
a state of quasi-particle and determining its energy accurate to the constant e n (0)iV. 



4 Green function for one-particle state 

It follows from Eq. (J23J) also that the quantity e n (ki) can be interpreted as an eigenvalue of the 
Hamiltonian for the quasi-particle excitation without taking into account interaction of quasi- 
particles. Therefore, the equation (J22J), written in the formalism of density matrix p^i. hk i = 
can be rewritten in the formalism of wave functions in coordinate representation and in 
a limit of large N, N — > oo in the following way: 

- (h+ + V sc )) Y.P [ Xr> = KmJ-e n (0))N6 rr ,, (25) 

' n 

where property J dk{n; kr\n; kr') = 5 rr > has been used; 5 rr > is a delta symbol. Since the 
energy e n (0) of the bound one-electron state is negative: e n (0) < 0, the right side of Eq. (J2SJ) 
represents itself a Dirac 5-function 8{r — r'). It allows to write Eq. ()25J) as: 

(4- hHF )T,&-, r ,r> = 6(r-r>), (26) 

where h HF = (h + if + V sc ), if = —J2 X . Eq. (|2l)|) is the equation for the Green function. It 
means that in the secondary quantized representation an operator 



possesses properties of non-perturbed Green function. 
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So, the quasi-particle excitation determined by the Hamiltonian h HF can be considered 
as a free particle whose equation of motion is the equation (j2HJ). 

In the many-body problem, in particular, in calculations of an energy-band crystal 
structure a contribution given by interaction of electromagnetic field with matter is played 
the essential role. To take into account many-particle effects due to a correlated motion 
of electron we should describe the system by self-consistent solutions of the non-stationary 
equation 



(2f 



where H is a Schrodinger hamiltonian in a non-relativistic case or a Dirac hamiltonian in 
a relativistic case. It turns the variational principle (jlUj) into a variational principle for 
the excited states, not being strictly definite one. Further, we shall show that within the 
framework of the concept of quasiparticle excitations this ambiguity can be removed by 
means of the account of interaction as, at first, a change of a quasi-particle mass and, 
second, as a changing of the location of reference point of energy band. 

We have proved that for the secondary quantized representation the operator p looks 
as p = and possesses properties of the Green function G\. Therefore the sum G\ over 

n from elements of the matrix p^, for the secondary quantized density matrix p describing 
an interacting particle satisfies a Dyson equation in a nonrelativistic case or to a Schwinger 
- Dyson equation in a relativistic case: 



Gi(l;2) = GS 0) (1;2)+ / d3 dA G^(l; 3)E(3, 4)Gi(4; 2 



(29) 



G[ (1;2) is a free Green function, £(3,4) is a self-energy operator: £ = XT + S c , E c is 
a correlation interactions, representing itself a part of the self-energy which describes the 
many-particle effects. Here numerical labels for the arguments are used: {ri,ti} — x± = 1, 
etc. Acting on Eq. (j2*9~)) by the operator ij^ — h HF and using the equation of motion for the 
free particle (}2l)j) we get the equation for the perturbed Green function as 



G x {n'\ 1, 2) - J rf3S c (n'; 1, 3)G 1 (n'; 3, 2) = 

= (-e n (0))N5 rir2 . 
Rewriting Eq. (|30p in the formalism of wave functions one gets 



(30) 



%— /; 

dt 



HF , 



ipnihri)- I c/f 2 S c (n; l,2)^ n (fcir 2 ) = 

(-e n (0))^ n (fciri). 



(31) 



Since the expression: i^j^- = e n (ki) takes place, then Eq. (|3*Tj) yields the Hartree - Fock 
taking into account of interacting quasi-particles 



h HF (ri)Mhrx) + I dr 2 E c (n;l,2)ij n (k ir2 ) 

(e n (0) + e n {k 1 ))^ n {k x r 1 ). 



(32) 
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Let us define a mass operator AM as: 

AMip n {k iri )= I dr 2 t c {n-l,2)^ n {k x r 2 ). (33) 



Within the framework of the concept quasiparticle excitations it is possible to represent 
the operator AM in the diagonal form: 

AM^hn) = (AM n (0) + AMnik^nikin). (34) 

It is known that the eigenvalue of mass operator possesses the property: AM n (ki) = 
AM n (—ki). Here AM n (0) is an eigenvalue of mass operator AM in the limit k — ► 0. 

It follows from here the physical sense of AM. It determines an effective mass of the 
quasi-particle and an efficient reference point of the energy band: 

h HF (r 1 )^) n (k l r 1 ) = (e n (0) + e n (&i))V>n(&i?"i) = 
= [(e n (0) + AM n (0)) + (e n (fci) + AM n (h))] Mhn). (35) 

Since the change of the mass of quasi-particle determined by the operator AM maintains 
the condition © then if to take into account the change of the reference point of band at 
interaction, the variational principle for the interacting system becomes a strictly definite 
one and takes the form: 

6E n $ n ] = 0, e n (0) = —ExtoE n (ki) - AM n (0). (36) 

Equation of motion for one-particle state and basis set of wave functions for exited atom 

Let us consider a Green function normalized per unit volume V = 1 so that an average 
energy in V is equal to the energy of the one-particle state and N = 1. If — e n (0) — > oo then 
Eq. (|3(J|) describes a propagation of one particle and should be rewritten as 



din'; 1, 2) - J d3t c (n'; 1, 3)Gi(n'; 3, 2) = 

= (-en(0))<W (37) 



It follows from here that according to the definition of Green functions we have the following 
expression for the energy e„(0): 

- e n (0) = C-a n , C -> oo; (38) 

where a n is a finite quantity Hence, since the energy is counted off from an arbitrary value, 
Eq. (}3l)j) yields the following expression for reference points e(0)^ of a quasi-particle energy 
and a antiquasi-particle energy 

e(0)± = ±a n = (ExtiEih) t AM n (0)) /2. (39) 

Here one took into account that N = 1; an extremum of zone is redefined as ExtrE(ki) = 
ExtrE n (ki) —C n , the sign {±} in left-hand side denotes a case of quasiparticles and antiquasi- 
particles, respectively; and the energy of particles in the pair is counted off from zero level. 
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One gets from the expression (J39|) that a n is the energy which is required to create a pair 
from quasiparticle and antiquasiparticle when k% = because 

a n = (e(0)+ - e(0)-)/2. (40) 

Because of an additional term e n (0) in a right-hand side of Eq. (|33|) we, generally 
speaking, cannot examine the left-hand side as a Hamiltonian operator of the quasi particles 
system acting on a corresponding wave function and as a consequence, can not construct 
a basis set of one-particle states of the problem. However, further we show, that h HF is a 
Hamiltonian of an electron - hole pair. 

Non-relativistic case 

One can examine in non-relativistic limit quantum systems which are characterized by 
a small value of AM n (0): 

AM n (0) -> 0. (41) 

It means that weak many-particle effects occur and, accordingly, we can speak about a 
"light" electron. The equality occurs under condition of (jJTJonly in the case if a n = 0. 
From here it follows, that the energy a n of a pair is equal to zero. In other words, the energy 
is not expended to create an electron - hole pair. 

Substituting Eqs. (J3~%j) . (}4~Tj) into Eq. and taking into account the condition a n = 0, 
one gets the Schrodinger equation as 

h HF ( ri )i> n (n) = tnMn), (42) 

which describes the quasiparticle - antiquasiparticle pair (a non-relativistic electron - hole 
pair). Here i n = e n — C. 

Since the energy a n , expended on creation of a pair, equals to zero we have proved 
that the variable e n can be understood as the energy of an electron - hole pair. Therefore, 
Eq. (jUJ) has a group of dynamic symmetry, which algebra is so(3) so(3) ~so (4) if to neglect 
an exchange interaction. As is known, a nonrelativistic hydrogen-like atom possesses such 
symmetry. Hence, we have proved that to calculate quasiparticle states in the non-relativistic 
case it is possible to use a basis set of states of a nonrelativistic hydrogen-like atom. 

However, for a heavy electron AM n (0) > 1 according to the formula (|40|) we always 

have 

a n = -AM n (0)/2 (43) 

and, hence, there does not exist equation such as Schrodinger one for its describing. From 
here we conclude that the heavy electron can not be examined in a nonrelativistic limit. 

Relativistic case 

Let us generalize the proposed approach to relativistic case. To do it we substitute 
Eqs. (13~%|) and (J4~3~j) into (|33j) and let n tends to n — ► oo: 

h^irJMkm) = (^MM. + ~ in (k x )^ Vfefari), n - oo. (44) 

Then, one can assume that the operator -§^ — h HF in Eq. (|44|) is a quasirelativistic hamiltonian 
written in the implicit form in the Hartee - Fock approximation. 

From consideration carried out above it follows that the desired relativistic equation of 
motion should describe a charged composite system from a pair of particles and have the 
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dynamic symmetry S0(4). A spin of given quantum system should be equal 1 as motion 
of a hole is a motion of an electron in many-particle positively charged matrix. In ^HJ the 
equation of motion of a relativistic charged vector boson has been found and shown, that 
it describes a relativistic hydrogen-like atom. The relativistic charged vector-boson appears 
a composite system with a corresponding spectrum of masses and in quasirelativistic limit 
n — > oo its energy E\ is determined by the expression: 



m 


rwy 2 


TTry 4 , 


' 4 




m7 6 / 


~2 


2n 2 


8n 3 \ 




n j 


8n 4 \ 



n 2 n\k\ k 2 



(45) 



Comparison of right-hand sides of formulas (|44J) and (|45)l yields that AMoo = lim^oo AM n (0) 
m is a rest mass m of an electron. Hence, Eq. (|4*4"j) is an equation of motion for a relativistic 
electron - hole pair with a reduced mass — m/2 which, apparently, is the relativistic 

charged vector-boson considered in quasirelativistic limit n — > oo. 



5 Conclusion 

So, the variation method to find interacting quasiparticle states in crystals was developed. 
The quasiparticle propagator was constructed by summation over elements of secondary 
quantized density matrix. This approach allows us to find motion equations of one-electron 
states of excited atom in crystals. 

References 

[1] D.R. Hartree. //Proc. Cambr. Phil. Soc. Vol.24, p.89 (1928). 
[2] V.A. Fock. //Zs. Phys. Bd. 61, s.126 (1930). 

[3] D. Hartree. Calculations of atomic sructures. (Foreign Literature, Moscow, 1960) 
[4] C. Froese-Fisher. The Hartree - Fock method for atoms. (Wiley, N.Y., 1977). 
[5] W. Kohn and L.J. Sham. Phys. Rev. A. Vol. 140, 1133 (1965). 

[6] Aryasetiawan F. and Gunnarson O. The GW method. / /Rep.Prog.Phys. 61. P.237-312 
(1998). 

[7] M.G. Veselov, L.N. Labsovsky. Theory of atom: Electronic shell structure. (Science, 
Moscow, 1986). 

[8] L.D. Faddeev, O.A. Yakubovsky. Lectures on quantum mechanics. (Leningrad univercity 
publishers, Leningrad, 1980) 

[9] P.O. Ldwdin. //Phys. Rev. Vol. 97, 1474 

[10] H.V. Grushevskaya, L.I. Gurskii. A projection operator technique for solution of 
relativistic wave equation on non-compact group: the case of a charged vector- 
boson. //Report BGUIR. Vol.1, no. 2, p. 12-20 (2003). E-print archive: www.arXiv.org, 



quant-ph/0301176 (2003) 



12 



[11] Slater J.C. Wave functions in a periodic potential. //Phys. Rev. (1937). Vol. 51, No. 
10. P.846-851 

[12] Wigner E., Seitz F. On the constitution of metallic sodium. //Phys. Rev. (1933). Vol. 
43, No. 5. P.804-810 

[13] H.V. Grushevskaya, L.I. Gurskii. Quantum-mechanical methods for calculation of elec- 
trophysical propeties of nanostructured systems. / /Report BGUIR. No. 2, p. 173 - 185 
(2004). 



